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Abstract—Optical coherence tomography (OCT) as retina 
imaging technology is currently used by ophthalmologist as 
a non-invasive and non-contact method for diagnosis of age- 
related degeneration (AMD) and diabetic macular edema (DME) 
diseases. Fluid regions in OCT images reveal the main signs of 
AMD and DME. In this paper, an efficient and fast clustering 
in neutrosophic (NS) domain referred as neutrosophic C-means 
is adapted for fluid segmentation. For this task, a NCM cost 
function in NS domain is adapted for fluid segmentation and then 
optimized by gradient descend methods which leads to binary 
segmentation of OCT Bscans to fluid and tissue regions. The 
proposed method is evaluated in OCT datasets of subjects with 
DME abnormalities. Results showed that the proposed method 
outperforms existing fluid segmentation methods by 6% in dice 
coefficient and sensitivity criteria. 

Index Terms—Optical coherence tomography; Neutrosophic 
theory; Fluid segmentation; Retina diseases 


I. INTRODUCTION 


Ophthalmologists use optical coherence tomography (OCT) 
as a non-contact method to diagnosis and follow up of 
retina disease. Fluid regions are detectable in OCT images 
and reveal abnormalities of retina disease including diabetic 
macula edema (DME) and age-related macular degeneration 
(AMD) (1). Image processing and data mining applications 
have been widely used in many application [2|-[7]. Fluid 
regions can be segmented automatically by image processing 
and data mining tools. Many fluid detection methods have 
been proposed in the literature such as kernel regression 
based segmentation of OCT images with DME (8). fluid-filled 
region delineation of boundaries in OCT (9). semi-automatic 
segmentation method for retinal OCT images tested in patients 
with diabetic macular edema (10), label propagation and 
higher-order constraint-based segmentation of fluid regions in 
OCT images of DME subjects (11). computerized assessment 
of intraretinal and subretinal fluid regions in OCT images 
of the retina (12). automated segmentation of pathological 
cavities in OCT scans [13], three-dimensional analysis of 
retinal layer and fluid-filled regions in OCT of the macula 
[14], probability constrained graph-search-graph-cut method 
for three-dimensional segmentation of fluid-associated abnor- 
malities (15), fluid segmentation with shortest-path methods 
and neutrosophic (NS) theory [16J-21], fluid segmentation 
with deep convolutional networks (CNNs) [22]-[25]. 
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Smarandache proposed NS theory in 1995, a branch of phi- 
losophy which studies the origin, nature and scope of neutrali- 
ties, as well as their interactions with different ideational spec- 
tra 26]. This theory has been applied on many applications 
including image segmentation (27-30). image thresholding 
|31], image edge detection (32), content-based image retrieval 
—(37], retinal image analysis [16}—(20], speech processing 
38], data clustering (29). and uncertainty handling (41), 
42]. The main contribution of this work is to optimize a cost 
function in NS domain to binarize retina OCT scans to fluid 
and tissue regions. The cost function is derived from fuzzy c- 
means (FCM) and is optimized by gradient descend methods 
in subsequence iterations. The rest of this paper is organized 
as follows: Section II reviews FCM and NS theory. Section 
III presents the proposed method. Experimental results and 
conclusion are discussed in sections IV and V, respectively. 


II. REVIEW OF NS AND FCM 


Since the cost function in this work is an extension of FCM 
and is defined in NS domain, here, a review on FCM and NS 
is discussed shortly. 


A. FCM 


The FCM is a clustering method that assigns a membership 
degree in interval [0-1] to each data point, Therefore, each 
data point is assigned to all clusters with different membership 
degrees and the sum of memberships to all clusters should be 
1. FCM cost function is defined as follows: 
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where u;; represents the membership degree of data point 
X; to the i” cluster with center vj. ||-|| and m are a norm 
metric and a constant variable for determining the fuzziness 
of the resulting partition, respectively. Minimizing FCM cost 
function leads to the following equations for the computation 
of membership degrees and cluster centers: 
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Firstly, membership degrees are initialized randomly and 
cluster centers are computed by Eq. (3). Then, membership 
degrees and cluster centers are computed repeatedly, until there 


is no significant changes of these parameters in subsequent 
iterations. 
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B. Neutrosophic 


NS theory is a newly brought up branches of philosophy. In 
this theory, set A, Anti-A (contemplates A in a correlation with 
its opposite), Neut-A (neutrality of A, neither A nor Anti-A) 
are defined (26). In each application, NS sets should be defined 
for all data points as well as relations between sets. For ex- 
ample, in our problem of interest as image segmentation, each 
pixel in NS domain is modeled as P(t, 7, f) meaning that it is 
t percent true, 7 percent indeterministic, and f percent false. 
Therefore, three sets Pys(i,7) = {T (4,7), 1(4, 9), Fi, 7) } are 
defined for pixels in NS domain as follows (27): 
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where, g(i,j) represents the intensity value of the pixel 
P(i,j), G is image matrix g filtered by average filter with 
window size w. Difference of two matrixes g and g is 
computed and stored in 6. 


III. PROPOSED METHOD 


In this work, NCM cost function is used and adapted for 
OCT fluid segmentation. Since OCT scans contain noise and 
there is ambiguity between fluid pixels and many pixels in 
background and tissue, indeterminacy set J in NS domain can 
be very useful and models such ambiguities easily. Here, each 
pixel in OCT scans are transferred to three sets T, J and F' by 
Eqs. (4), (6) and (8), respectively. Then, NS sets are presented 
for clustering cost function as follows: 
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Gradient descent method is used for cost function optimiza- 
tion which leads to the following relations for the computation 
of T’, J and F’ sets in NS domain and cluster centers. 
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Parameter AK is a common part in Eqs. (13)-(I5p. It is 
computed once and used three times which leads to speedup in 
each iteration of optimization. In the first step, 12 clusters are 
considered for clustering. The reason is that there is 11 layers 
in retina and each layer can be assigned to one cluster based on 
its texture and gray level. One extra cluster is considered for 
fluid regions as fluid cluster. After convergence, clusters are 
sorted ascendingly based on the gray level of cluster centers. 
The first cluster is considered as fluid cluster with label | and 
other clusters are considered as tissue clusters with label 0. 
Therefore, the proposed clustering scheme leads to a binary 
segmentation of OCT scans. 


IV. EXPERIMENTAL RESULTS 
A. Evaluation Metrics 


To show the effectiveness of the proposed method, it has 
been evaluated with three metrics; dice coefficient, precision 
and sensitivity; computed from true positive (TP) fluid pixels 
detected as fluid, false positive (FP) tissue pixels detected as 
fluid, true negative (TN) tissue pixels detected as tissue and 
false negative (FN) fluid pixels detected as tissue. 
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These criteria are used to evaluate fluid segmentation re- 
sults of automated methods in comparison with manually 
segmented fluid regions by ophthalmologist experts. 


B. Dataset 


The proposed fluid segmentation method is tested on a 
dataset from OPTIMA Cyst Segmentation Challenge which 
contains 4 subjects with 49 images per subject where the image 
resolution varies from 512 x 496 to 512 x 1024. The fluid 
regions of each OCT image are manually segmented by two 
ophthalmologist experts as ground truth images. This dataset 
is publicly available and can be found onlind' 


C. Resultes 


Segmented fluid regions by the proposed method are com- 
pared with ground truth images. Fig. [1] depicts fluid regions 
segmented by expert and the proposed method for two scans, 
each scan in one row. The proposed fluid segmentation method 
is compared with 4 fluid segmentation methods proposed in 
and [44-46]. The results of these methods are shown in 
Fig. |2| It is clear visually that the proposed method has lower 
number of false negaives which leads to higher accuracy. 


(a) (b) (c) 


Fig. 1: Fluid segmentation results for 2 scans: (a) Input OCT 
image, (b) manually segmentated by expert and (c) 
segmented by the proposed method. 


Table [I|reports the dice coefficient, sensitivity and precision 
of the proposed method and methods in and [44]-(46]. 
Reported results in Table |I/ are the comparison between all 
methods and fluid regions manually segmented by expert 1. 
Table [II|reports the same results but in comparison with expert 
2. In comparison with the results of expert 1, the proposed 
method outperforms other methods with respect to dice co- 
efficient and sensitivity of 82.23% and 90.59%, respectively. 
These measures of the proposed method in comparison with 
expert 2 are 81.62% and 90.11% which also outperforms 
other methods. For precision measure, methods in with 
85.29% and with 90.87% have the best performance in 
comparison with expert | and 2, respectively. 


V. CONCLUSION 


This paper presented a fluid segmentation method based 
on NCM cost function in NS domain. Minimizing the cost 
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Fig. 2: Segmented fluid by: (a) expert, (b) method in 17}, 
(c) method in [44], (d) method in and (e) the proposed 
method 


TABLE I: Dice coefficient, sensitivity and precision of the 
proposed method and other methods in comparison with expert 
1. 


Expert 1 Sub Method | Method | Method | Method | Prop. 
" ) in in in in Method 

1 73.49 80.4 71.40 82.96 83.17 

73.90 55.10 45.49 78.11 74.61 

Dice Coeff | 3 78.46 75.35 69.54 82.23 84.85 
4 78.12 71.78 TLAS 80.75 86.32 

Ave 75.99 70.66 64.39 81.01 82.23 

1 70.81 82.19 72.49 84.43 86.95 

2 96.79 99.04 70.45 98.94 98.42 

Sensitivity 3 13.72, 85.13 47.38 85.18 86.58 
4 78.78 80.59 T1.719 84.49 90.51 

Ave 80.52 86.73 67.02 88.26 90.59 

1 93.00 85.06 54.87 84.03 82.12 

74.36 54.18 31,12 78.48 69.31 

Precision 3 94.89 79.88 30.93 85.45 83.43 
4 96.97 88.62 54.98 93.20 75.01 

Ave 89.80 76.93 47.97 85.29 TTA6 


function resulted in binary segmentation of Oct images into 
tissue and fluid regions. Segmentation results on Optima 
dataset showed that the proposed method outperforms other 
segmentation methods in dice coefficient and sensitivity mea- 
sures while for precision criteria, other methods had the best 
performance. Future efforts will be directed towards proposing 
a cost function in NS domain for better modeling of noise and 
uncertainty in OCT pixels. NS can model uncertainty of OCT 
pixels in deep convolutional networks and leads to more robust 
network against noise and imaging device. Therefore, using 
NS theory for fluid segmentation by CNN can be proposed as 
another future work. 
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